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ABSTRACT 

Dissipation of magnetohydrodynamic (MHD) wave energy has been proposed as a viable heating mechanism 
in the solar chromospheric plasma. Here, we use a simplified one-dimensional model of the chromosphere to 
theoretically investigate the physical processes and the spatial scales that are required for the efficient dissipa¬ 
tion of Alfven waves and slow magnetoacoustic waves. We consider the governing equations for a partially 
ionized hydrogen-helium plasma in the single-fluid MHD approximation and include realistic wave damping 
mechanisms that may operate in the chromosphere, namely Ohmic and ambipolar magnetic diffusion, viscosity, 
thermal conduction, and radiative losses. We perform an analytic local study in the limit of small amplitudes to 
approximately derive the lengthscales for critical damping and efficient dissipation of MHD wave energy. We 
find that the critical dissipation lengthscale for Alfven waves depends strongly on the magnetic field strength 
and ranges from 10 m to 1 km for realistic field strengths. The damping of Alfven waves is dominated by 
Ohmic diffusion for weak magnetic field and low heights in the chromosphere, and by ambipolar diffusion 
for strong magnetic field and medium/large heights in the chromosphere. Conversely, the damping of slow 
magnetoacoustic waves is less efficient, and spatial scales shorter than 10 m are required for critical damping. 
Thermal conduction and viscosity govern the damping of slow magnetoacoustic waves and play an equally 
important role at all heights. These results indicate that the spatial scales at which strong wave heating may 
work in the chromosphere are currently unresolved by observations. 

Subject headings: Sun: atmosphere — Sun: chromosphere — Sun: magnetic fields — Sun: oscillations — 
waves 


1. INTRODUCTION 

The plasma heating of the upper solar atmosphere is one 
of the long-standing problems in solar physics. The physical 
processes responsible for the transport of energy from the so¬ 
lar interior and its dissipation in the atmospheric pl asma are 
under intense research (see, e.g. JParnell & De Moortell2012l) . 
One of the mechanisms that has been proposed to explain the 
transport and dissipation of energy involves the propagation 
and damping of ma gnetohyd rodynamic (MHD) waves (see 
the recent review bv lArreguH|2015i and references therein). 
Indeed, observations indicate that MHD waves are ubiquitous 
in the solar atmosphere and can have the required ene rgy to 
heat the upper layers (see, e.g ., iCargill & d e Moortel J20J1J; 
iMcIntosh et al.ll20TTt lHahn & Savinll2014t Hess et all 20151) . 
While the overwhelming presence of the waves is demon¬ 
strated by the observations, the physics behind the damping of 
the waves and the deposition of wave energy into the plasma 
remains poorly known. 

In this paper, we perform a theoretical study on the condi¬ 
tions for efficient dissipation of MHD wave energy in the solar 
chromosphere. The realistic modelling of the physical pro¬ 
cesses in the chromospheric plasma is challenging. Because 
of the relatively low temperature, the chromospheric plasma 
is only partially ionized and neutrals are predominant at low 
and medium altitudes. It has been shown that partial ioniza¬ 
tion effects have a strong impact on chromospheric dynam¬ 
ics (see, e.g JMartmez-Svkora et al. 1120121 [Leake et al.ll2014l) . 
Ion-neutral collisions may play a crucial role in the release of 
magn etic energy in the form of heat (iKhomenko & Colladosl 
120121) . Therefore, the consideration of partial ionization is a 
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unavoidable requisite for the realistic description of the chro¬ 
mospheric physics. 

The role of partial ionization on the damping of Alfven 
waves has been investigated in detail in the literature. Es¬ 
timations of the damping rate due to ion-neutral collisions 
as a function of height in the chromosphere indicate that 
the damping is most efficient for waves with high fre¬ 
quencies closer to the local ion-neutral collision frequency 
(see, e.g ., |De Pontieu et all 120011: [ khodachenko et al.l l2004t 


ILeake et all 120051 ISoler et all I2012al among others). Re- 
centlv. ISoleretalJ (120151) investigated in detail this phe¬ 
nomenon and showed that high-frequency Alfven waves can 
be critically damped, i.e., overdamped, because of ion-neutral 
collisions. The energy carried by these overdamped waves 
could be efficiently deposited in the plasma as a result of 
the strong dissipation. In fact, previous computations of the 
plasma heating rate obtained from numerical simulations (see 
Goodmanll201 U ISong & Vasvliunasll201 It iTu & Son~gll201 3t 


iRussell & Fie tcherl 1201 3l) showed that dissipation of Alfven 
wave energy can provide a sustained heating over time that 
is sufficient to compensate the chromospheric radiative losses 
at low altitudes. So, it is well established that partial ioniza¬ 
tion effects are important for the correct study of Alfven wave 
damping and associated heating in the chromosphere. 

Concernin g magnetoacoustic waves, Khoda chenko et afl 
(120041120061) obtained that thermal conduction and viscos¬ 
ity c an be impor t ant in their dam ping. Als o, results from, 
e.g., iPorter et all (119941) and iCarbonell et al.l (120041) suggest 
that the effect of radiative losses should be taken into ac¬ 
count depending on the plasma physical conditions. Con¬ 
versely, ion-neutral collisions may play a relatively less im¬ 
portant role in the direct damping of magnetoacoustic waves 
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(see iForteza et al .1120071) . However, partial ionization should 
be appropriately accounted for, since the presence of neutrals 
in the plasma modifies the efficiency of thermal conduction 
and viscosity. Therefore, as discussed by Kho dachenko et ai] 
(j2 006i ), the correct and complete description of the damping 
of the various types of MHD waves in the solar atmosphere 
in general, and in the chromosphere in particular, requires the 
consideration of all energy dissipation mechanisms, includ¬ 
ing thermal, collisional, and frictional effects. We add that the 
correct and complete description of the damping also requires 
the consideration of all the components of the plasma, namely 
the various ionized and neutral species. 

Here we are concerned about the physical mechanisms that 
are actually relevant in damping the waves and about the 
spatial scales that are required for the chromospheric damp¬ 
ing to be efficient. We introduce the concept of ‘critical 
dissipation lengthscale’, which is related to the occurrence 
of wave cutoffs. The existence of cutoffs caused by the 
strong damping of the waves is well documented in the litera- 
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ers). The physical nature of the wave cutoffs resides in the fact 
that the damping is so strong that the wave restoring force is 
effectively suppressed (see a more extensive explanation in 
ISoler et aPl2013bl) . As a result, wave perturbations decay in a 
timescale much shorter than the wave period, a phenomenon 
called critical damping or overdamping. For practical pur¬ 
poses, this means that the waves are unable to propagate and 
so they cannot transport energy away from the chromospheric 
medium. All the energy initially stored in the wave pertur¬ 
bations is eventually deposited in the plasma in the form of 
heat. Therefore, overdamped waves are strong candidates to 
be connected with efficient plasma heating in the chromo¬ 
sphere ISoler et al.ll2015l) . Wave cutoffs due to critical damp¬ 
ing are physically and conceptually different from the cutoffs 
due to plasma stratification studied in many previous works. 

The purpose of this paper is twofold. On the one hand, we 
aim to perform a comprehensive theoretical study of the phys¬ 
ical mechanisms capable of producing the cutoff of Alfven 
and magnetoacoustic waves under chromospheric conditions. 
It is common among theoretical works that investigate wave 
damping to focus on the influence of one specific damping 
mechanism and ignore other effects. We are interested in pro¬ 
viding a consistent description of all relevant physical pro¬ 
cesses that may be involved in strong wave damping. On the 
other hand, we aim to obtain the critical dissipation length- 
scales for the various types of MHD waves. The value of 
the critical dissipation lengthscale is relevant because it deter¬ 
mines the spatial scales needed for wave heating to become 
efficient. Strong plasma heating produced by dissipation of 
MHD wave energy would necessarily require spatial scales of 
the order of the critical dissipation lengthscale. 

Another important question from the theoretical point of 
view is the role of helium in the damping of the waves. The 
majority of previous works that investigated wave damping 
ignored the influence of helium and assumed a chromospheric 
plasma composed of hydrogen alone. However, it has been 
shown that the presence of neutral and ionized helium can en¬ 
hance the damping of Alfven waves due to ion-neutral colli¬ 
sions (IZaaarashvili et al.1l2013l) . A secondary objective of this 
work is to determine the impact of helium on the efficiency 
of the damping mechanisms and, consequently, on the critical 
dissipation lengthscales. 


This paper is organized as follows. Section [2] contains the 
description of the chromospheric model, the wave dissipation 
mechanisms, and the basic equations governing wave propa¬ 
gation and damping. The critical dissipation lengthscales for 
Alfven waves and slow magnetoacoustic waves in the chro¬ 
mospheric model are investigated in Sections[2and|4] respec¬ 
tively. Later, the results of this article are discussed in Sec¬ 
tion [5] and some concluding remarks are given in Section [6] 
Finally, the effect of dissipation on the nonlinear coupling 
of Alfven waves and slow waves is briefly explored in Ap¬ 
pendix [a] 


2. BASIC EQUATIONS 
2.1. Partially ionized chromospheric model 

We adopt a one-dimensional, static, gravitationally strati¬ 
fied model for the chromosphere based on the semi-empirical 
model F of iFontenla et al'l (119931) . hereafter the FAL93-F 
model. The reason for choosing the FAL93-F model in¬ 
stead of more recent models is that the models tabulated in 
IFontenla et al.l (1993) explicitly provide the variation of the 
neutral and ionized helium densities with height. The same 
model was used in the previous work by ISoler et afl (120151) . 
but the influence of helium was not considered there. 

We treat the chromospheric medium as a partially ion¬ 
ized hydrogen-helium plasma composed of electrons, pro¬ 
tons, neutral hydrogen, neutral helium, and singly ionized he¬ 
lium. Hereafter, subscripts e, p, H, He I, and He II explicitly 
denote these species, respectively. The presence of doubly 
ionized helium. He III, and that of heavier species is ignored 
because of their negligible abundance in the chromosphere. 
We define the fraction of species p = e, p, H, He I, or He II as 



(1) 


where pp = nipnp is the mass density of species p , with nip 
and up the mass particle and number density, and p = 2/; Pp 
is the total mass density. We note that p e « 0 owing to the 
very small electron mass. We assume a strong thermal cou¬ 
pling and use the same temperature, T, for all the species. Fig- 
ure[Tfa)-(c) shows the dependence on height, /?, of the temper¬ 
ature, total density, and fraction of the various species, respec¬ 
tively. The sharp chromosphere-to-corona transition region is 
located at h ~ 2,000 km, where the temperature increases and 
the density decreases abruptly. Hydrogen is mostly neutral for 
It < 1,500 km, becomes to be ionized for h > 1,500 km and 
is fully ionized for h > 2,000 km. We note that there is a rel¬ 
atively large abundance of neutral helium at large altitudes in 
the chromosphere because the temperature is not high enough 
to fully ionize helium. 

The various species in the plasma exchange momentum by 
means of particle collisions. The strength of the interaction 
depends on the so-called friction coefficient, a. The friction 
coefficient for collisions between two charged species q = e, 
p, or He II and q' - e, p, or He II, is (e.g.. ISpitzerl 1 1 962b 
Bra ginskiill 19651) 
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where m qq - = m q m q > / (m q + m q ) is the reduced mass, e is the 
electron charge, is Boltzmann’s constant, eo is the permit¬ 
tivity of free space, and In A qq f is Coulomb’s logarithm given 













































Wave heating scales in the chromosphere 


3 


by (e.g. JSnitzerlll 962HVranies & Krstidl2013l) 


In A m 


= In 
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In turn, the friction coefficient for collisions between a 
charged or neutral species /? = e, p, H, He I, or He II, and 
a neutral speci es n = H or He I is (e.g.. iBraginskiil IT965b 
IChapman & Cowlingll 1970b 


&/3n — 


8 k B T 


nmp n 


1/2 


Cpu, 


(4) 


where crg n is the collision cross section. Recent estima tions of 
this parameter can be found in lVranies & Krsticl (120131 ) and an 
extensive discussion on the importance of its value for wave 
damping is given in lSoler et al . (201 5). 

The expressions of the friction coefficients given above are 
also valid for self-collisions, i.e., collisions between particles 
of the same species. The total friction coefficient of species /3 
with the other species is 

Qp = 'y, oipp'i (5) 

P'*P 


while the total friction coefficient of species f3 including self¬ 
collisions is ap t tot = <*p + app. 


2.2. Single-fluid MHD description of plasma dynamics 

We study the dynamics of the partially ionized chromo¬ 
spheric plasma within the framework of the single-fluid MHD 
approximation. The single-fluid MHD approximation as¬ 
sumes a strong coupling between ions, electrons, and neu¬ 
trals so that all the species effectively behave as one fluid (see 
iBraginskiiil 1965b . This approximation is valid in the chro¬ 
mosphere as long as the wave frequencies remain lower than 
the ion-neutral collision frequency of the plasma and the 
wavelengths remain longer than the mean free path of par¬ 
ticles between collisions (see, e.g., Zaaara shvili et al.l 1201 lL 
120131 ISoler et al.ll2013bl) . In this approximation, the basic 
MHD equations are written in terms of average quantities, 
while the effects of the interactions, i.e., collisions, between 
the various spec ies remain in the form of nonideal terms (see 
details in, e.g., [ ktrteza et all 120071 iMeier & Shumlakl 120121 
Khomenko et al. 1201 4) . The basic single-fluid MHD equa¬ 
tions used in this work are 
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( 8 ) 
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( 10 ) 


with p the mean atomic weight of a hydrogen-helium plasma 
given by 


1 
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P ~ 2£p + (TH + T^Hel + T^HeH 
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Equations (l 6 l) (ITOt are the continuity equation, momentum 
equation, induction equation, energy equation, and the equa¬ 
tion of state, respectively. In these equations, ^ = Jy + v • V 
denotes the material or total derivative, p is the thermal pres¬ 
sure, v is the velocity vector, B is the magnetic field vec¬ 
tor, g is the acceleration of gravity, y is the adiabatic index, 
p the magnetic permeability, n is the viscosity tensor, rj and 
/ 7 a are the coefficients of Ohmic and ambipolar diffusion, re¬ 
spectively, £. represents the net effect of all the sources and 
sinks of energy, and R is the ideal gas constant. We note that 
Hall’s term is ignored in the induction equation due to its ir- 
relevant role in d amping the waves (Zaaaras hvili et al.1l2012l 
ISoler et alJUoOb . 


2.3. Dissipation mechanisms 

Here we discuss the importance of the various dissipation 
mechanisms that appear in the governing Equations (l 6 l) (flOt . 

Equation ([ 8 ]> contains two magnetic diffusion terms. The 
first term on the right-hand side of Equation ([Sj is Ohmic dif¬ 
fusion, which is caused by electron collisions. The coefficient 
of Ohmic diffusion is given by 


a e 

1 2 ’ 
A ie l n% 


( 12 ) 


Equation (IT 2 l > corresponds to the classic Spitzer’s coefficient 
of magnetic diffusion ( ISpitzerll 1962b . The second term on the 
right-hand side of Equation © is ambipolar diffusion. This 
term contains the effect of ion-neutral collisions. The coeffi¬ 
cient of ambipo lar diffusion in a hydr ogen-helium plasma is 
given by (see lZaaarashvili et al.ll2013b 


I + + 2^H^He I^H He I 

7a = --7-j-7-• d3) 

I ]-j e j J 

Figure[l}d) displays the ratio |B | 2 i]a/p as a function of height 
in the chromospheric model for three different values of the 
magnetic field strength. We see that ambipolar diffusion is 
dominant throughout the chromosphere except at low heights, 
where Ohmic diffusion can be of importance for weak mag¬ 
netic fields. Therefore, we shall consider both Ohmic and 
ambipolar diffusion in our computations. 

Viscosity in a partially ionized plasma is essentially deter¬ 
mined by self-collisions of ions and neutrals. Electron vis¬ 
cosity can be s afely neglected by virtue of the small elec tron 
mass (see, e.g., IBraginskiil965biMeier & ShumlaU2012b . On 
the one hand, we assume the neutral viscosity tensor, 7 r n , to 
be isotropic because of the absence of the direct effect of 
the magnetic field on neutrals. Isotropy of neutral viscosity 
might be altered by ion-neutral coll isions (seelVranies| |2014 ), 
but this effect is usually neglected (IMeier & Shumlak|l2012i) . 
Thus, the neutral viscosity tensor used here is 


= ~fn | Vv + (Vv) T - /|V • V 


(14) 


where / is the identity tensor and is the coefficient of neutral 
viscosity given by 




'W H 7I 2 
^ ^H,tot 


'»HeI»h e r 

^He I,tot j 


(15) 


where both neutral hydrogen and neutral helium are included. 
On the other hand, the ion viscosity tensor, iq, has a compli¬ 
cated form in the presence of a magnetic field. It is usually 
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Figure 1 . Variation of physical parameters with height above the photosphere according to the chromospheric model F of lFontenla et ali j 1993 '): (a) temperature, 
(b) total density, (c) fraction of species, (d) ratio of ambipolar to Ohmic diffusivities for three different magnetic field strengths, (e) ratio of neutral viscosity to 
ion parallel viscosity, and (f) ratio of neutral thermal conductivity to electron parallel thermal conductivity. 


described as the sum of five components accounting for par¬ 
allel (or compressive) viscosity, perpendicular (or shear) vis¬ 
cosity, and gyroviscosity (see the full expression in lBraginskiil 
119651) . In the magnetized chromospheric plasma parallel vis¬ 
cosity is several orders of magnitude larger than both perpen¬ 
dicular viscosity and gyroviscosity. Then, we only consider 
the ion parallel viscosity component, namely 




-3fy 




: Vv, 


(16) 


where b = B/ |B| is the unit vector in the direction of the mag¬ 
netic field and is the coefficient of ion parallel viscosity 
given by 


&II = 0.96 


'«p«p fWHeIl»HeII^ 


0(p,tOt 


#He II,tot 


knT, 
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where both protons and singly ionized helium are included. 
Figure|T]e) displays the ratio £n/<Ti,|| as a function of height in 
the chromospheric model. We clearly see that neutral viscos¬ 
ity is several orders of magnitude more important than ion 
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viscosity throughout the chromosphere. Because of this re¬ 
sult, we are allowed to completely neglect ion viscosity com¬ 
pared to neutral viscosity. We shall approximate ft « n n in our 
computations. 

The effect of all the sources and sinks of energy is included 
in the function X in the right-hand side of the energy equation 
(Equation (0). The expression of X is 

X = -V • q - L (p, T) + j • E* - ;r: Vv + //, (18) 


where the various terms on the right-hand side are: the di¬ 
vergence of heat flux due to thermal conduction q = -kVT, 
with k the thermal conductivity tensor; the radiative loss 
function L(p, 7); the generalized Joule heating j • E*. with 
j = (V x B) Ip the current density and E* = E + v x B the 
effective electric field; the viscous heating n : Vv; and an ad¬ 
ditional and unspecified source of heating 77. 

As happens for viscosity, thermal conduction in a magne¬ 
tized plasma is strongly anisotropic. For convenience, we split 
the thermal conductivity tensor into its components parallel, 
Aj|, and perpendicular, k ± , to the magnetic field direction. In a 
fully ionized plasma k\\ is essentially determined by electrons, 
while k± is mainly due to ions and is negligible. In a partially 
ionized plasma the isotropic thermal conductivity of neutrals 
has to be added to both parallel and perpendicular compo¬ 
nents. Thus, the parallel and perpendicular components of 
thermal conductivity are approximated by k\\ » /c c || + a:,, and 
k ± ~ A- n , where k c .\\ and K n are the parallel electron thermal 
conductivity and the total neutral thermal conductivity, re¬ 
spectively, given by 


, = 3.2 


nlkjT 

^e,tot 


Kn 3 


, O'H^ot OfHe,tot j 


k\T. 


(19) 

( 20 ) 


FigureHJe) displays the ratio y as a function of height in 
the chromospheric model. The neutral thermal conductivity 
is found to be more important than the electron thermal con¬ 
ductivity except at large altitudes, where both thermal con¬ 
ductivities are of the same order of magnitude. This is so be¬ 
cause hydrogen is largely ionized at those large altitudes, but 
helium remains mostly neutral. Hence, we shall retain both 
neutral and electron thermal conductivities in the computa¬ 
tions and express the divergence of the heat flux as 


V ■ q = — B • V . 

|B| 2 


b ■ vr - v ■ (/r n vr). 


( 21 ) 


The radiative loss function, L ( p , T), accounts for plasma 
cooling owing to radiative losses. Determining the chromo¬ 
spheric radiative losses as function of temperature and den¬ 
sity is a difficult task that requires complicated numerical so¬ 
lutions of the radiative transfer problem. The radiative loss 
rate depends, e.g., on the completeness of the atomic model 
used for the calculation, on the atomic processes included, on 
the ionization equilibrium, and element abundance assumed, 
among other factors. The full solution of the radiative trans¬ 
fer problem is beyond the purpose and scope of the present 
paper. Here we follow a frequent alternative approach to ac¬ 
count for the solar plasma radiative losses based on a semi- 
empi rical parametrization of the radiative loss function (see. 




l2001tlSchure et al 


iRosner et al.1 Il978t iKlimchuk & C argill 


20091) . This method enables us to incorpo¬ 


rate radiative losses in an approximate way without the need 


of solving the full radiative transfer problem. The inconve¬ 
nience of this approach is that the semi-empirical radiative 
loss function is obtained under the assumption of optically 
thin plasma, while the cool chromospheric plasmas of interest 
here do not completely satisfy this condition. Owing to fi¬ 
nite optical thickness, the actual radiative losses of the plasma 
would be probably reduced in some degree compared to the 
optically thin proxy. This fact should be taken into account 
when interpreting the effect of radiative losses on the results. 
The expression of the radiative loss function we use is 


L(p, T)=p 2 x*T a , (22) 

whereof* and a are piecewise constants depending of the tem¬ 
perature. We use the parametrization of x* and a given in Ta¬ 
ble 1 of ISoler et al . ( 2012b ). which are obtained from up-to- 
date computations of radi ative losses derived from the CHI- 
ANTI atomic database dDere et all 1 19971 lLandi et afl 1201 2l) 
assuming typical ab undances in the sola r at mosphere. The 
reader is referred to iParenti et alJ (120061) and iParenti & Viall 
(2007) for details. 

The generalized Joule term j • E* takes into account plasma 
heating because of dissipation of electric currents. The ex¬ 
pression of j ■ E* is 

j - E* = yt/77 |j N | 2 -t- yU77c |jj_| 2 , (23) 

where 77 c is the so-called Cowling’s diffusivity given by 

77 c = 77 + |B| 2 77 A, (24) 

and j|| and j, are the components of the current density paral¬ 
lel and perpendicular to the magnetic field, respectively, that 
can be computed as 


j|| = — [(V x B) • b] b, 

(25) 

p 

j_L = -b x [(V x B) x b]. 

(26) 


B 


Thus, Ohmic magnetic diffusion is caused by the dissi¬ 
pation of parallel currents, while Cowling’s magnetic dif¬ 
fusion, i.e., the joint effect of Ohmic and ambipolar diffu¬ 
sion, is cau sed by the dissipation of perpendicular currents 
(iKhomenko & Colladosll2012l) . As explained before, we shall 
retain both Ohmic and ambipolar diffusion in the following 
computations. 

2.4. Equations for one-dimensional wave propagation 

We reduce the set of basic equations to the case of the one¬ 
dimensional chromospheric model. We use a Cartesian co¬ 
ordinate system and assume that the plasma properties vary 
along the x-direction only, whereas y and z are ignorable di¬ 
rections. Therefore, the x-direction corresponds to the ver¬ 
tical direction, so that gravity is oriented in the negative in¬ 
direction, namely g = (-g, 0,0), with g the acceleration of 
gravity at the solar surface. We can conveniently rotate the co¬ 
ordinate system for v and B to lie on the xy-plane with no loss 
of generality. We redefine the x and y directions as the paral¬ 
lel, ||, and perpendicular, ±, directions, respectively. Hence, 
v = (v||,Vj_,0) and B = (By, 7J ± , 0). The solenoidal condi¬ 
tion V • B = 0 imposes that dB\\/dx = 0, while from the x- 
component of Equation 0 we get dB\\/dt = 0. Therefore, By 
is a constant in both space and time. In this one-dimensional 
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case, Equations I© © become 
dp d 

a = "5t (p,|l) ' 
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(27) 

(28) 

(29) 

(30) 




+ H 


(31) 


In turn, the variations of temperature are related to those of 
pressure and density by 

L— = (32) 

T dt p dt p dt 

Let us consider the static case so that we set vy = v ± = B ± — 
0 and d/dt — 0. From Equation ( |28| > we get 


dp 

dx 


= ~Pg- 


(33) 


This is the condition for a gravitationally stratified plasma and 
is satisfied by the chromospheric model. In turn, from Equa¬ 
tion lOTT ) we find that the condition for the plasma to be in 
thermal equilibrium is that the arbitrary heating function must 
be 


H — L(p, T) - 


d_ 

dx 


, ,dT 

^ + Kn) dx 


(34) 


The arbitrary heating term essentially balances the back¬ 
ground radiative losses since the contribution of thermal con¬ 
duction, i.e., the second term on the right-hand side of Equa¬ 
tion (l34l >. is almost negligible. We note that this heating term 
is merely included to formally maintain the semi-empirical 
chromospheric model in thermal equilibrium. The presence 
of this arbitrary heating have no effect on the behavior and 
damping of the waves superimposed on the background chro¬ 
mosphere. 

Inspection of Equations (|27]>—(f3Tb reveals that the velocity 
longitudinal component, vy, and the thermal pressure, p, are 
nonlinearly coupled to the velocity and magnetic field trans¬ 
verse components, v ± and B ± . This means that Alfven waves, 
which are polarized in the transverse direction to the magnetic 
field, can nonlinearly drive perturbations associated with lon¬ 
gitudinally polarized slow magnetoacoustic waves. On the 
contrary, if Vj_ and B ± are initially zero, slow magnetoacous¬ 
tic waves cannot drive Alfven waves in the plasma. 


2.5. Approximate local analysis of perturbations 

The purpose of this paper is to obtain expressions for the 
lengthscales that govern the dissipation of MHD waves in the 


chromospheric plasma. To do so, we consider local pertur¬ 
bations superimposed on the background plasma and perform 
an approximate study in the limit of small amplitudes. The 
use of the local analysis is justified by the fact that, for the 
waves to be dissipated in the chromosphere, the wavelengths 
and the associated dissipation lengthscales must be necessar¬ 
ily shorter than the chromospheric gravitational scale height. 
The pressure scale height in the chromosphere is ~ 300 km, 
so that we must consider shorter lengthscales in this analy¬ 
sis. Consequently, the effect of gravity on the wave perturba¬ 
tions is ignored in the present local analysis, although gravi¬ 
tational stratification is fully taken into account in the back¬ 
ground plasma. 

We define the dimensionless parameter e = max \B ± \ /By 
and assume that e <s 1. Then we write 


P ~Po + eV, 

(35) 

P-Po + e 2 p'. 

(36) 

r = r 0 + e 2 r. 

(37) 

v ll = e2v ll, 

(38) 

Vj. = ev' ± , 

(39) 

B\\=B 0 , 

(40) 

B ± = eB' ± , 

(41) 


where the subscript 0 denotes a background quantity and the 
prime ' denotes a perturbation. To separate the perturbations 
of Alfven waves from those of slow waves, the perturbations 
of the perpendicular components of velocity and magnetic 
field are assumed to be first-order in e, while the perturba¬ 
tions of the remaining quantities are assumed to be second 
order in e. We substitute these quantities in Equations (l27l >- 
(|3TT > and separate the various terms according to their order 
with respect to e. 

The conditions for static equilibrium (Equations (l33l ) and 
GS) are consistently recovered from the zeroth-order terms 
in e. The first-order equations in e govern the behavior of v x 
and B ± , and so they describe linear Alfven waves, namely 


dv' ± - d 2 v' ± | B 0 dB' ± 
dt n dx 2 ppo dx 


dB' ± dv 1 

■ =D i 


d 2 B\ 


dt 


0 - 5 — + t]c , , 
dx dx~ 


(42) 

(43) 


where = f B /po- The second-order equations in e involve 
v'|, //, //, and r, so that they describe linear slow magne¬ 
toacoustic waves and also the generation of slow magnetoa¬ 
coustic waves due to the nonlinear coupling with the Alfven 
waves, namely 


dp' _ ovy 

dv \\ _ dp' , 4<r n d 2 Vy 1 dB' 2 

P ° dt dx 3 dx 2 2 p dx 


(44) 

(45) 


dp' 

~dt 


= -ypa 


d Z i 

dx 


(7-1) 


d 2 T' 

tail + *n) —ffT ~ P0 L pP' 


, 1 IdB'V (dv' 

-PoLtT' + -W:bH +&hr 
/i \ OX I \ ox 


(46) 


where L p and L r are the partial derivatives of the radiative loss 
function with respect to density and temperature, respectively. 
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namely 


L p - 


(f) , 

\ d P) PO ,To 


L T = 



(47) 


In addition, p', p', and T are related through the equation of 
state as 



Po Po T 0 


(48) 


The equations for higher orders in e represent nonlinear cor¬ 
rections on the Alfven and slow waves. In this approximate 
study we consider sufficiently low amplitudes for the high- 
order corrections in e to be negligible. Therefore we restrict 
ourselves to the first-order and second-order equations in e. 


3. CRITICAL DISSIPATION LENGTHSCALE OF 
ALFVEN WAVES 


To start with, we consider Alfven waves, which are gov¬ 
erned by Equations (142} and (|43} . These two equations in¬ 
volve the transverse components of v and B and can be ap¬ 
propriately combined into an equation for B' ± only, namely 


d 2 B’ ± 2 d 2 B' ± 
dt 2 A dx 2 


'+<fn 


d 3 B\ 


+ PcZ 


d 4 B\ 


wher e c\ = B 2 /ppo is the Alfven velocity squared. Equa¬ 
tion (1491 ) governs linear Alfven waves damped by Cowling’s 
diffusion and viscosity. It can be shown that v' ± satisfies the 
same equation. 

Let us assume that a packet of Alfven waves is generated 
impulsively by a perturbation with a certain spatial form. Any 
disturbance in the plasma can be expressed as a superposition 
of Fourier modes. Each Fourier mode is characterized by a 
wavenumber, 4, and a frequency, a>(k). Thus, we write 


B' ± = ^ exp [ikx - i<o{k)t ], (50) 

k 


where Bi- is the amplitude of the Fourier mode with wavenum¬ 
ber k. We substitute this expression into Equation ((49} to 
find the dispersion relation for small-amplitude, linear Alfven 
waves in a viscous-resistive plasma, namely 


cj 2 (k) + ik 2 (t/ c + fn) cj(k) -k 2 c\- k A p c ln =0. (51) 

Equati on (l5~il w as previously derive d by, e.g. JChandrasekharl 
(1196 111 and Zaaaras hvili et all (1 201 2l) . 

The individual Fourier modes can be studied separately 
since no interaction between different Fourier modes of the 
wavepacket takes place in the small-amplitude, linear regime 
considered here. We focus on one individual Fourier mode 
of the wavepacket, so we set k to a fixed value and oj(k) is 
obtained by solving Equation (BH . The analytic solution to 
Equation (dTT) is 


Lo(k) = ±kc\ 


1 - 


k 2 (PC ~ <Tn) 

4c 2 


2 \ 


1/2 


- iff (pc + lr) ■ (52) 


The first term on the right-hand side of Equation (152} is the 
real part of the frequency, namely <air(4), and the second term 
is the imaginary part, namely a>i(k). The + and - signs in 
front of wr (k) correspond to waves propagating towards the 
positive A-direction (upwards) and towards the negative x- 
direction (downwards), respectively. Hence, the initial per¬ 
turbation naturally splits into two wavepackets propagating 


in opposite directions. The imaginary part of co(k) defines a 
damping time scale for the Fourier mode, so that its ampli¬ 
tude decreases in time following an exponential law, namely 
exp (- |wi(4)| t). Ideal undamped Alfven waves are recov¬ 
ered when no dissipation is present ( 77 c - l n - 0), so that 
wr (k) = ±kcA and a>i(k) = 0 in that case. However, in 
the presence of Cowling’s diffusion and/or viscosity Alfven 
waves are damped. 

The quality factor, Q(k), is a dimensionless parameter that 
characterizes how efficiently damped a specific Fourier mode 
is. The quality factor is commonly defined as 




wr (k) 

a>i(k) 


(53) 


The quality factor is independent of the direction of wave 
propagation. The strength of wave damping depends on the 
value of Q(k). According to the definition of Q(k) given in 
Equation d53} . the waves are weakly damped or underdamped 
when Q(k) > 1/2. In that case, most of the wave energy 
would not dissipate in the chromosphere. The larger Q(k), 
the less efficient damping, so that no wave energy dissipa¬ 
tion takes place if Q(k) —> 00 . Conversely, when Q(k) <1/2 
the waves are overdamped and the wave energy dissipation is 
very strong: most of the wave energy is dissipated. The most 
extreme situation takes place when Q(k) = 0, which corre¬ 
sponds to a wave cutoff. In a cutoff scenario, waves cannot 
propagate and all the energy stored in the waves is deposited 
in the plasma. The expression of Q(k) computed from the 
Alfven wave frequency (Equation (152} ) is 


Q(k) = 


k-4n| -ffk 1 

r h 

0 


if k < kA. 
if k > kA, 


(54) 


where kA is a cutoff wavenumber for Alfven waves given by 


k& — 


2c a 


17C - <fn| 


(55) 


The Fourier modes of the wavepacket with k > k A cannot 
propagate in the form of Alfven waves since their quality fac¬ 
tor is zero. As discussed before, the physical reason for this 
result is that damping is extremely efficient for these large-4 
modes, so that their propagation is inhibited due to the very 
strong energy dissipation. Thus, the wavepackets that would 
eventually propagate away from the chromosphere necessar¬ 
ily contain Fourier modes with 4 < 4 a alone, while the Fourier 
modes with 4 > 4 a would be completely dissipated in the 
chromospheric plasma. 

The cutoff wavenumber defines the lengthscale for which 
dissipation is total, namely 


4a - 


2 n nVrjc-lA 


kA 


ca 


(56) 


The critical dissipation lengthscale is plotted in Figure (2ja) 
as a function of height for three values of the magnetic field 
strength, namely Bo = 10, 100, and 1000 G. In the lower chro¬ 
mosphere, the larger the magnetic field strength, the shorter 
the critical dissipation lengthscale. The behavior is the oppo¬ 
site one in the middle and upper chromosphere, where the crit¬ 
ical dissipation lengthscale increases when the magnetic field 
strength is increased. The reason for this result resides in the 
mechanisms responsible for the damping of the waves. This is 
explored in Figure |2}b), which shows the critical dissipation 
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(a) (b) 




(4) 




Fi gure 2. (a) Critica l dissipation lengthscale for Alfven waves as a function of height above the solar photosphere according to the chromospheric model F of 
IFontenla et alj <19930 . The various linestyles correspond to a magnetic field strength, Bo, of 10 G (solid), 100 G (dotted), and 1000 G (dashed), (b) Effect of the 
various damping mechanisms on the critical dissipation lengthscale for Alfven waves with Bo = 10 G: full result (solid), effect of Ohmic diffusion alone (dotted), 
effect of ambipolar diffusion alone (dashed), and effect of viscosity alone (dot-dashed), (c) Threshold magnetic field strength for which ambipolar diffusion 
dominates over Ohmic diffusion as a function of height, (d) Effect of helium on the critical dissipation lengthscale for Alfven waves with Bo = 10 G: full result 
(solid) and result without helium (dotted). 


lengthscale for Bq = 10 G computed when only one specific 
damping mechanism is retained. We recall that Cowling’s dif- 
fusivity, rjc, contains both Ohmic, //, and ambipolar, //,\, dif- 
fusivities. It can be seen that Ohmic diffusion dominates the 
damping at low heights, while ambipolar diffusion becomes 
the most important mechanism at large heights. Viscosity 
plays no important role throughout the chromosphere. The 
switch from Ohmic damping to ambipolar damping depends 
on the value of the magnetic held strength. This is shown 
in Figure|2jc), where the threshold magnetic held strength is 
plotted as a function of height. This threshold magnetic held 
strength is computed as B 0 ~ ^jr]/ij A . We see that strong 
helds of the order of kG are required at the low chromosphere 
for ambipolar diffusion to be more important than Ohmic dif¬ 
fusion. However, only a few tens of G are needed in the 
medium chromosphere, and very weak helds are needed in 
the upper chromosphere. These results concerning the relative 
importance of Ohmic and ambipolar diffusion are consistent 
with the expected importance of the two effects according to 
Figure [ljd). 

We have also explored the particular effect of helium on 
the value of the critical dissipation lengthscale. Figure [2fd) 
compares the Alfven wave critical dissipation lengthscale for 
Bo = 10 G when helium is taken into account and when the 


effect of helium is neglected in the dissipation mechanisms. 
This hgure indicates that the effect of helium is not impor¬ 
tant when the damping is caused by Ohmic diffusion. For 
the magnetic field strength used in Figure |2fd), this happens 
when h < 900 km. However, helium should be taken into ac¬ 
count when the damping is dominated by ambipolar diffusion. 
In Figure |2}d), this happens when h > 900 km. The reason 
for the important impact of helium on the ambipolar damping 
is that helium remains largely neutral in the chromosphere. 
Obviously, since the importance of helium is linked to that 
of ambipolar diffusion, the particular height at which the ef¬ 
fect of helium becomes important depends on the value of the 
magnetic held strength considered (see again Figure [2}c)). In 
general, the influence of helium should not be neglected (see 
also lZaqarashvili et al.ll2013l) . 


4. CRITICAL DISSIPATION LENGTHSCALE OF SLOW 
MAGNETOACOUSTIC WAVES 

We move to the slow magnetoacoustic waves and consider 
Equations (l44t — (f46l) . They can be combined to arrive at a 
rather involved equation for v' alone, namely 


ar iy - 1 Ta?- W 


d_ 

dt 


^i^L\ d A 

3 dx 2 ) dt 
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(a) (b) 




Figure 3. ( a) Critical dissi pation lengthscale for slow magnetoacoustic waves as a function of height above the solar photosphere according to the chromospheric 
model F of IFontenla et ali <19930 . The various linestyles correspond to t he fu ll result (solid), the effect of thermal conduction alone (dotted), the effect of 
viscosity alone (dashed), and the analytic approximation given in Equation ED (symbols), (b) Effect of helium on the critical dissipation lengthscale for slow 
magnetoacoustic waves: full result (solid) and result without helium (dotted). 




dx 2 dt 


+ (y- 1) 




dx 2 




dt 


- (y - 1) U 


dx 2 


~ CJj 


|2 d '2 


d l B 


dxdt 


— (y — l) 


dxdt 


— 1 

(d_B\\ 

PPO 

l dx ) 


dv ^ 
dx 


(57) 


where c 2 = 'ypo/po is the sound velocity squared and k , cj p , 
and cot are defined as 


- 

K= -(K e ,|| +*n), 

(58) 

P 0 


9 


- P °T 

OJp — E/p, 

(59) 

po 


PoT 0 T 

OJj =- L*t • 

(60) 


P o 


The general solution to Equation (l57l > is the sum of the so¬ 
lution of the related homogeneous equation and the partic¬ 
ular solution of the inhomogeneous equation. On the one 
hand, the solutions of the homogeneous equation are uncou¬ 
pled slow magnetoacoustic waves. The uncoupled slow waves 
do not interact with the Alfven waves studied in Section [3] 
On the other hand, the particular solution of the inhomoge¬ 
neous equation represents the generation of slow magnetoa¬ 
coustic waves due to the nonlinear coupling with the Alfven 
waves. We note that the right-hand side of Equation (l57l > in¬ 
volves terms with B' ± and v' ± . These terms are associated with 
the Alfven waves and are responsible for nonlinear driving 
slow waves. Here, we focus on the solution of the homo¬ 
geneous version of Equation (l57l > that represents uncoupled 
slow waves. The nonlinear driving of slow waves due to the 
coupling with the Alfven waves is briefly addressed later in 
Appendix IaI 

We perform a Fourier analysis and write vjj as 

v[| = ^ < V\\,k exp [ikx - iu(k)t ], (61) 

k 


where 'Vp is the amplitude of the k -th Fourier mode of v' 
and k and oj(k) denote again the wavenumber and angular fre¬ 


quency. We obtain the dispersion relation from the homoge¬ 
neous version of Equation (l57l >. namely 

a?(k ) + ici 2 tJ 2 (k) - aicj(k) — iao — 0, (62) 

with 

02- + (y~ l)(kk 2 + l>t), (63) 

a\-k 2 [c 2 + (y - \)(kk 2 + co T J , (64) 

a 0 = (y - 1) | kk 2 + y (oj t - w p )j k 2 . (65) 

Equation (l62l) is a third-order polynomial and is the dispersion 
relation of parallel-propagating slow magnet oacoustic waves, 
but it also describes thermal modes dFieldl I1965h . In gen¬ 
eral, slow waves and thermal modes are coupled when nona- 
diabatic effects are c onsidered (see, e.g.. iPorter et al.lfl994l : 
iCarbonell et alj|2004l) . 

The thermal mode is the nonadiabatic version of the so- 
called entropy mode, which is related to nonpropagating per¬ 
turbations of density in a plasma (Goedbloed & Poedts 2004). 
Following the same method as in iSoler et aTI (!2012bb , an ap¬ 
proximate solution for the thermal mode in the case of weak 
damping can be found by neglecting terms of O (w 2 j in Equa¬ 
tion ( l62l >. By doing so, we find 


(y- 1 )[kk 2 + j (oj T - OJ p ) 

a>(k) ~ -i ---r- 

c 2 + (y- 1 )(kk 2 + oj t ) 


( 66 ) 


The approximate thermal mode frequency given in Equa¬ 
tion d66l > is purely imaginary. This means that thermal modes 
are always nonpropagating regardless of the spatial scale of 
their perturbations. For chromospheric conditions, the ther¬ 
mal mode imaginary part of the frequency is always negative, 
so that_the thermal mode is a purely damped solution (see, 
e.g.. ISoler et all[2012bl> . Growing modes related to thermal 
instabilities would require a positive imaginary part of the fre¬ 
quency. Therefore, according to Equation (l66l) thermal insta¬ 
bility does not take place in the chromospheric plasma. Ther¬ 
mal instability requires higher temperatures and is though to 
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occur in the coronal plasma (see details in iFieldlTl965!) . For 
an optically thin plasma, the condition for thermal instabil¬ 
ity is that a < 0 , where a is the exponent of the temperature 
in Equation (l22t . In the absence of thermal conduction and 
radiative losses. Equation (f66l) gives oj(k) = 0. This zero- 
frequency solution is the classic entropy mode in an adiabatic 
plasma (see lGoedbloed & Poedtsll2004) . We shall not explore 
thermal modes further in this paper. 

We go back to the study of slow waves. The exact an¬ 
alytic solution to Equation (l62t is too complicated to ob¬ 
tain a simple expression of the critical dissipation length- 
scale of slow waves. However, we can apply the method of 
[Soler et all (l2013blfah and use the polynomial discriminant of 
Equation (l62l) to obtain some useful information. To do so, we 
perform the change of variable a>(k) = is(k) in Equation (l62t 
and remove the common factor -i. Then a cubic equation in 
s with real coefficients is obtained, namely 


damped because of the effect of radiative losses alone. The 
impact of radiative losses on the damping of slow waves in the 
chromosphere is therefore negligible. Conversely, the critical 
dissipation lengthscale is very similar to the full result when 
either thermal conduction or viscosity are considered alone 
(these results are overplotted in Figure[3ja)). The two mecha¬ 
nisms are equally important and both should be considered to 
explain the slow wave dissipation. 

Taking advantage of the numerical result displayed in Fig¬ 
ure | 2 a), we can try to find an analytic approximation of the 
slow wave critical dissipation lengthscale. The process is as 
follows. First, we decouple the slow wave from the thermal 
mode by roughly neglecting the independent term from Equa¬ 
tion (l62l i. Then, we consider the numerical result that radia¬ 
tive losses are unimportant, so that we are allowed to neglect 
the terms with <x>j and co p . This gives an approximate solution 
for the slow wave frequency as 


s(k ) 3 + a. 2 s(k ) 2 + ciis(k) + gq - 0 (67) 

Subsequently, we compute the polynomial discriminant. A, of 
the cubic equation using the standard formula as 

A = a\a\ - 4a\ - 4a\ao - 21 + l 8 a 2 aiciQ. ( 68 ) 


a>(k) =3 ±kc s 


4c? 


4£. 

3 


(y- l)k 


1/2 


k 2 (4{ 0 


+ (y- l)k , 


(69) 


The full developed expression of the discriminant is not given 
here for the sake of simplicity. The discriminant has the gen¬ 
eral property that it is zero when the original polynomial has 
a multiple root. This is very convenient, since a multiple zero 
of Equation d62l) occurs precisely when the slow wave suf¬ 
fers a_cutoff (see an explanation for the reason of this result in 
ISoler et ak|[2013blal) . As a function of height we numerically 
determine the real values of k, i.e., the cutoff wavenumbers, 
which cause the discriminant to vanish. We find that the dis¬ 
criminant has only one real root, namely k s , which defines the 
slow wave critical dissipation lengthscale as ,l v = 2n/k s . The 
numerically obtained critical dissipation lengthscale of slow 
waves is plotted in Figure [3j a) as a function of height. The 
lengthscale varies several orders of magnitude from ~ 10 4 m 
in the lower chromosphere to ~ 10 m in the upper chromo¬ 
sphere. Thus, the critical dissipation lengthscale for slow 
waves is significantly shorter than that obtained for Alfven 
waves in Section [3] Slow waves require very short spatial 
scales to be efficiently damped. Another difference between 
Alfven waves and slow waves is that the slow wave critical 
dissipation lengthscale is unaffected by the value of the mag¬ 
netic field strength. 

We note that the slow wave critical dissipation lengthscale 
reaches very short values in the lower chromosphere. These 
extremely short lengthscales are probably stretching the as¬ 
sumptions behind the theory used in this paper. Such short 
lengthscales are of the same order of or shorter than the 
mean free path of particles between collisions. Therefore, 
the single-fluid MHD approach may be compromised for such 
small scales in the lower chromosphere. Multi-fluid or kinetic 
approaches would be more convenient instead. Readers must 
be aware that the results obtained here for the slow wave crit¬ 
ical dissipation lengthscale in the lower part of the chromo¬ 
sphere should be interpreted with caution. 

In order to gain physical insight on the mechanisms respon¬ 
sible for the slow wave cutoff, now we compute the criti¬ 
cal dissipation lengthscale when only one specific damping 
mechanism is retained. In the case of radiative losses, we 
find that this mechanism cannot produce the cutoff of the slow 
wave. In other words, the slow wave never becomes critically 


where again the + and - signs in front of the real part of the 
frequency correspond to upward and downward propagating 
waves, respectively. We use Equation (l69l ) to compute quality 
factor, Q(k), namely 

( |%-(y-l)*| -Jkr-k 1 

(?0) 

0 if k>k s , 


where k s is the approximate cutoff wavenumber for slow 
waves given by 


2c s 


R 1 - (y- l)k 


(71) 


The critical dissipation lengthscale defined by this cutoff 
wavenumber is 


As 


2 n ~ 

ks c s 


(72) 


We overplot Equation (f72l> in Figure[3ja). An excellent agree¬ 
ment between the approximate lengthscale and the numeri¬ 
cally computed lengthscale is obtained. In addition, we can 
compare Equation (l72l > with the expression for the Alfven 
wave critical dissipation lengthscale (Equation [56l>. We see 
that both equations are formally similar and that thermal con¬ 
duction plays in slow waves the same role as Cowling’s diffu¬ 
sion plays in Alfven waves. 

As done in the case of Alfven waves, we have also inves¬ 
tigated the specific effect of helium on the critical dissipa¬ 
tion lengthscale for slow magnetoacoustic waves. Figure [3th) 
shows the slow wave critical dissipation lengthscale when he¬ 
lium is taken into account and when the effect of helium is 
omitted from the dissipation mechanisms. We obtain that the 
influence of helium is negligible except at large heights. As 
already discussed, the reason for this result is that helium re¬ 
mains mostly neutral at all heights in the chromosphere, while 
hydrogen gets fully ionized at large heights. Therefore, at 
those large heights there is a significant amount of neutral he- 
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Figure 4. Critical dissipation lengthscale of Alfven waves (solid line) as a 
function of height above the sola r pho tosphere for a magnetic field strength 
that varies according to Equation (73) with B p h = 2 kG. The critical dissipa¬ 
tion lengthscale of slow waves (dashed line) is also plotted for comparison. 

lium that efficiently contributes to both thermal conduction 
and viscosity. 

5. DISCUSSION 

Strongly damped waves are good candidates to produce sig¬ 
nificant heating of the chromospheric plasma via conversion 
of wave energy into thermal energy. Thus, it is important to 
know the physical processes involved in the damping and the 
spatial scales at which these processes work. We have in¬ 
vestigated the physical mechanisms that may be able to effi¬ 
ciently dissipate MHD wave energy in the chromosphere and 
have computed the spatial scales associated to critical damp¬ 
ing (wave cutoffs). 

By comparing the results of Alfven waves (Figure [2} and 
slow magnetoacoustic waves (Figure [3j we notice several im¬ 
portant implications. For instance, the physical mechanisms 
responsible for the occurrence of the cutoffs are different for 
the two types of MHD waves. The Alfven wave cutoff owes 
its existence to the effect of Cowling’s (Ohmic+ambipolar) 
diffusion, while viscosity plays a much less important role. 
For practical purposes the impact of viscosity on Alfven 
waves can be safely neglected. In addition, neither thermal 
conduction nor radiative losses affect the Alfven wave damp¬ 
ing at all, owing to the incompressible nature of Alfven waves. 
Conversely, the slow magnetoacoustic wave cutoff is caused 
by the combined effects of viscosity and thermal conduction, 
with the two mechanisms playing an equally important role at 
all heights in the chromosphere. Since slow waves are unaf¬ 
fected by the magnetic field Cowling’s diffusion does not play 
a role in the cutoff of slow waves, and radiative losses are 
again irrelevant. Among all the considered damping mech¬ 
anisms, radiative losses is the only mechanism that can be 
safely neglected on the case of both Alfven and slow waves. 
These results are consistent with previous estimations of the 
efficien cy of various damping mechan isms in the chromo¬ 
sphere (iKhodachenko et aUl2004l.r2006l) . 

Another relevant result is that the values of the critical dis¬ 
sipation lengthscales are rather different for Alfven and slow 
magnetoacoustic waves. In the case of Alfven waves, we 
find that the critical dissipation lengthscale is very sensitive 
to the magnetic field strength (see Figure [2f a)). In the low 
chromosphere the stronger the field, the smaller the dissipa¬ 
tion lengthscale. Conversely, in the middle and high chromo¬ 


sphere the stronger the field, the larger the dissipation length- 
scale. Because of the simplicity of our ID model, we were 
restricted to consider straight and constant magnetic fields. 
In reality the chromospheric magnetic field is neither straight 
nor constant but composed of flux tubes that expand from the 
photosphere to the corona, with the magnetic field strength 
varying from strengths of kG in the photosphere to strengths 
of few G in the corona. To incorporate the effect of a magnetic 
field that realistically varies with height, let us assume that the 
chromospheric magnetic field strength can be roughly approx¬ 
imated by the semi-empirical formula used bv lLeake & Arberl 
(2006), namely 

/ \ 0 ' 3 

B = B p h (—) , (73) 

VPph/ 

where B p h and p p h are the magnetic field strength and den¬ 
sity at the photosphere. We compute the Alfven wave critical 
dissipation lengthscale according to this prescription for the 
magnetic field strength. This result is displayed in Figure [4] 
For this realistically varying field strength, the Alfven wave 
critical dissipation lengthscale is roughly constrained in be¬ 
tween 10 m (low chromosphere) and 1 km (medium and heigh 
chromosphere). We should note that the Alfven wave critical 
dissipation lengthscale obtained for a varying field strength is 
roughly of the same order of magnitude as that obtained by 
ISoler et ail (120151. see their Figure 3), although in that paper 
the effect of helium was ignored. Conversely, the slow magne¬ 
toacoustic wave critical lengthscale is unaffected by the mag¬ 
netic field strength. For comparison purposes, we have over¬ 
plotted in Figure[4]the slow wave critical lengthscale that was 
already shown in Figure[3ja). We see that the slow wave crit¬ 
ical lengthscale is several orders of magnitude shorter than 
that obtained for Alfven waves. This comparison informs 
us that the dissipation mechanisms working in the chromo¬ 
spheric plasma are more efficient in dissipating magnetic en¬ 
ergy (associated to Alfven waves) than acoustic energy (asso¬ 
ciated to slow magnetoaocustic waves), because the efficient 
dissipation of slow waves requires shorter spatial scales. 

As mentioned in the Introduction, the observational evi¬ 
dence of MHD waves in the chromosphere is overwhelming. 
However, the fact that waves are observed does not automati¬ 
cally imply that those waves are actually contributing in heat¬ 
ing the plasma. According to the results of this work, only 
the waves with the appropriate spatial scales (wavelengths) 
may be able to efficiently damp and deposit their energy in the 
chromospheric medium. These spatial scales, below 1 km for 
Alfven waves and even smaller for magnetoacoustic waves, 
are unresolved by current space-based and ground-based tele¬ 
scopes. In this direction, future instruments operating at very 
high temporal and spatial resolutions like, e.g., the Atacama 
Large Millimeter/submillimeter Array (ALMA), may be nec¬ 
essary to observe those very short spatial scales and so to 
understand how MHD waves _actually contribute to chromo¬ 
spheric plasma heating (IWedemever et al.B2015lf . 

6 . CONCLUDING REMARKS 

To conclude with, it is worth mentioning the possible im¬ 
pact of the simplifications used in this work. We have de¬ 
fined the critical dissipation lengthscales as the spatial scales 
at which strict wave cutoffs occur. A strict cutoff is math¬ 
ematically defined so that the wave quality factor vanishes, 
i.e., Q = 0. However, there are several physical effects not 
considered in the present investigation that could remove the 
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mathemat ically strict wave c utoffs. IZaaarashvili et al.1 (120121) . 
and later ISoler et aTI (120151) . showed that Hall’s term in the 
induction equation ha s the effect of rem oving the strict cut¬ 
offs. As explained in IZaqarashvili et al.ll (t2012h . in the pres¬ 
ence of Hall’s term the waves do not suffer a strict cutoff when 
the wavenumber exceeds the critical value. Instead, the qual¬ 
ity factor is nonzero but is always constrained in the interval 
0 < Q < 1/2. This is the interval corresponding to over¬ 
damped waves. In connection to the amount of wave energy 
that can be deposited in the plasma, the fact that the waves 
have mathematically strict cut-offs (Q - 0) or are overdamped 
(0 < Q < 1/2) makes no practical difference. In the case of 
overdamped waves, the imaginary part of the frequency (re¬ 
lated to damping) is larger than the real part of the frequency 
(related to propagation). Therefore, as happens in a strict cut¬ 
off scenario, an overdamped MHD wave is unable to propa¬ 
gate, and so it cannot transport its energy away from the chro¬ 
mosphere ( ISoler et al.ll2015l) . Hence, the actual condition for 
wave energy to be efficiently dissipated in the chromosphere 
should be Q < 1 /2, which is less restrictive than the condi¬ 
tion Q = 0 imposed here. In both cases, the critical length- 
scales are very similar. Another effect that can remove the 
strict cutoffs is the consideration of electron inertia_terms in 
the multi-fluid descript ion of the plasma (Zaaarashvi li et all 
1201 2t ISoler et alll2015l) . As in the case of Hall’s term, elec¬ 
tron inertia terms replace the strict cutoff by overdamping, but 
such a result has no practical implications concerning wave 
heating. 

We have assumed that chromospheric plasma dynamics can 
be studied under the single-fluid MHD approximation. Multi¬ 
fluid ef fects on the wave cu t offs ha ve been investigated in de¬ 
tail by ISoler et al.l ([201 3bllal 120151) . They showed that wave 
cutoffs consistently occur in both single-fluid and multi-fluid 
theories. However, in the multi-fluid description it is found 
that waves may be able to propagate again when the length- 
scale is further reduced and beco mes sm aller than a second 
critical value. As explained by ISoler et al.l (l2013bl) . this sec¬ 
ond critical lengthscale would correspond to the spatial scale 
at which neutrals decouple from ions, and so the single¬ 
fluid approximation breaks down. Therefore, the existence of 
the second critical lengthscale is not captured by the single¬ 
fluid description of the plasma. This second critical length- 
scale could be so small for chromospheric conditions that 
the fluid treatment of the plasma may become compromised 
(ISoler et al.l 120151) . Also, a more appropriate treatment of 
the slow wave critical dissipation in the lower chromosphere 
would require multi-fluid or kinetic approaches because of the 
extremely short spatial scales obtained. 

We have considered a static one-dimensional model of the 
chromosphere, which allows us to analytically investigate 
the damping of Alfven and slow MHD waves in the limit 
of small amplitudes and for spatial scales satisfying the lo¬ 
cal approximation. The assumption of local analysis is fully 
justified in view that the critical dissipation lengthscales of 
both Alfven and slow waves are much smaller than the pres¬ 
sure scale height in the chromosphere (~ 300 km). Nev¬ 
ertheless, it should be acknowledged that a static model is 
a dramatic simplification because it misses the dynamic be¬ 
havior of the chromosphere seen in high-resolution obser¬ 
vations and reproduced by advanced numerical simulations 
(e.g.. iMartmez-Svkora et alll2012l) . A time-dependent back¬ 
ground should be necessary to account for the time-varying 
dynamics of the actual chromosphere. In addition, because of 
the linear, small-amplitude regime studied here, our approx¬ 


imate analysis is able to describe wave damping but misses 
the changes in the background due to wave energy deposition 
in the plasma. In order to properly overcome the limitations 
of the present analytic study, this paper needs to be extended 
in the future by considering fully nonlinear numerical simu¬ 
lations of MHD wave excitation and dissipation, including a 
time-dependent dynamic background and the self-consistent 
changes in the plasma due to wave energy dissipation. This is 
an interesting task and will be tackled in a forthcoming work. 
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APPENDIX 

EFFECT OF DISSIPATION ON THE COUPLING OF ALFVEN AND SLOW WAVES 


Here we briefly investigate the slow waves that are nonlinearly driven because of the presence of Alfven waves. To do so, we 
must explore the solutions to the full Equation (157b , taking into account the inhomogeneous driving term on the right-hand side. 
We are mostly interested on the impact of the Alfven wave dissipation on the process of mode coupling. Therefore, we shall 
focus on this particular aspect of the coupling between the two waves. A detailed analysis of the full evolution of nonlinear waves 
is left for future numerical studies. 

The analytic study of the inhomogeneous Equation (l57l > is difficult due to the complexity of the equation. In order to make 
further progress and to study the wave coupling analytically, for simplicity we shall omit all dissipative mechanisms from Equa¬ 
tion (l57b . By doing so, we neglect the effect of dissipation on the nonlinearly driven slow waves. However, the effect of dissipation 
on the primary Alfven wave is still considered at full. Under this approximation. Equation (l57l) simplifies to 

2 <3 2 v[, i g 2 B' 2 

_!i -r 2 _- =_±. (Ait 

dt 2 s dx 2 2/jpo dxdt 


To study the nonlinear coupling between the primary Alfven wave and the nonlinearly driven slow wave, we write B’ ± in the 
form of Equation (l50b . In turn, we write v' as in Equation (|6H . but now we use K and Q to denote the wavenumber and frequency 
of slow waves. We use this different notation to distinguish the wavenumber and frequency of slow waves from t he w avenumber 
and frequency of Alfven waves, namely k and a>. First, we substitute the expressions of B'_ and vj’i into Equation (lAlb and obtain 
the relation between the wavenumbers and frequencies of the primary Alfven wave and the nonlinearly driven slow wave, namely 
K — 2k and Q = 2a>. Therefore, we recover the typical result that the frequency and wavenumber of the nonlinearly generated 
slow mode are twice the frequency and wavenumber of the primary Alfven mode. We shall use these relations to write all the 
following expressions in terms of k and io. Next, we find the relation between the amplitudes of B' ± and as 


•Viut = 


kcj(k) 


ar(k) - k 2 cl 


2 BP 


(A2) 


where u>(k) is the Alfven wave frequency given by Equation ( l52l ). 

It is more useful to write the relation between the amplitudes of two components of velocity, v' ± and Vy. To do so, first we must 
find the relation between v' ± and B' ± . Hence, we write 


4 = k exp [ikx - iu>(k)t + /(Oo(k)], 

k 


(A3) 


where r V± J k is the amplitude of the Fourier mode of v' ± with wavenumber k, and <po(k) accounts for phase differences between 
and B' ± . If A < k A the relation between A’/ and 'V ± k can be cast as 


while <po(k) is given by 


c A 

Tu = -± 8 k . 

Bo 


ifo(k) = n ± arctan 



(A4) 


(A5) 


where, as before, the + and - signs correspond to upward and downward propagating waves, respectively. Therefore, the relation 
between the amplitudes of v^_ and B' ± i s co nstan t, wh ile the phase shift depends on k. The phase shift is n for k <sc k,\ and +n/2 
when k —> k A . Now, we use Equations (lA2b and (lA4b to write the relation between the amplitudes of two components of velocity 
as 


•Viut = 


kat(k) 


aj 2 (k) - k 2 c 2 


<y 2 

k l ,k 


(A6) 


In general, the relation of amplitudes of the two velocity components is a function of k. When k /c A , the damping of Alfven 
waves is very weak and we can approximate the Alfven wave frequency by tu(k) * ±kc A . Then, the relation of amplitudes of the 
two velocity components simplifies to 


'V 



(A7) 
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(a) (b) 




Figure 5. (a) Mode conversion height in the chromosphere as a function of the magnetic field strength in the dissipationless limit (k <§; k\). (b) Amplitude of 
the longitudinal component of velocity, 'Vp, as a function of height in the chromosphere. Symbols denote results in the dissipationless limit (k k&) and lines 
are the results in the limit of total dissipation (k —> k\). Three different values of the magnetic field strength are considered: 10 G (solid line and crosses), 100 G 
(dotted line and diamonds), and 1000 G (dashed line and asterisks). In all cases, we used < V±,k = 1 km/s as reference. 


which is independent of k. Equation (IA7b shows that the parallel velocity component associated to the nonlinearly driven slow 
waves diverges when c 2 A = c 2 . This is the condition of mode conversion from Alfven waves to slow waves. Since ca is a function 
of the magnetic field strength, the specific height in the chromosphere at which the condition c A = c 2 is satisfied depends on 
the value of Bq. Figure [2 a) shows the mode conversion height as a function of the magnetic field strength. For weak fields, 
the condition c 2 A = c 2 happens in the medium and/or upper chromosphere, while the mode conversion height moves to the low 
chromosphere when strong magnetic fields are considered. We have chosen 'V± t k = 1 km/s as a reference value for the amplitude 
of the transverse component of velocity and have plotted in Figure |5jb) t he a mplitude of the longitudinal component for three 
different values of the magnetic field strength. As predicted by Equation (1A7I) . the longitudinal component of velocity diverges 
at those specific heights where c A = c 2 . It is at those heights that the nonlinear driving of slow waves is most effective. In the 
dissipationless limit ( k <s k\), energy from the Alfven waves may be efficiently transferred to the slow waves at those specific 
heights. 

Conversely, in the limit k —> kp, the Alfven waves are fully damped and their frequency approximates by u>(k) ~ 
- ik 2 (r]c + in) /2. In such a case, the relation of amplitudes of the two velocity components tends to 


•Vn 


k ~ r 


\ric-ln\ 


c A 


Uc-?, 


I)' 


-°v 


_L ,k- 


(A8) 


c\+cl 


We have used Equation (lA8b to overplot in Figure [5fb) the amplitude of the longitudinal component of velocity in the limit of 
strong diffusion of Alfven waves (k —> k\). We see that the infinite amplitudes obtained in the difussionless case are replaced by 
finite amplitudes. To estimate the maximum amplitude of we set c 2 A = cl and consider that the effect of viscosity is much 
less important than that of Cowling’s diffusion for the damping of Alfven waves. Hence we neglect viscosity from Equation (1A8I) . 
Then, the equation becomes 

maxCVp)-T^Vl*- (A9) 

4c a 


Equation (lA9b indicates that for small-amplitude Alfven waves with r V± t k ^ ca, the maximum amplitude of the longitudinal 
component of velocity is necessarily very small in the limit of strong dissipation. 

The approximate results obtained here suggest that the nonlinear coupling of Alfven and slow waves is affected by dissipation 
of the primary Alfven wave. Dissipation may reduce the efficiency of Alfven-to-slow mode conversion in the chromospheric 
layer where c 2 A = c^. However, these results are obtained under quite restrictive approximations and their validity should be 
checked by means of fully nonlinear numerical simulations. We leave this task for forthcoming works. 











